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The two-leg Hubbard ladder is characterized by the ratio 
of the inter- to intra-leg hopping t±/t, the relative interaction 
strength U/t and the electron filling. Here, using density ma- 
trix renormalization group and Monte Carlo simulations, we 
examine the dependence of the pairing correlations on these 
parameters. We find that the pairing correlations are en- 
hanced when the top of the bonding quasiparticle band and 
the bottom of the antibonding band are near the Fermi level. 
We present results on the single-particle spectral weight and 
the antiferromagnetic correlations in order to explain this be- 
havior. 



One goal of numerical many-body calculations has 
been to explore the qualitative properties of a given 
model, in particular, the nature of the dominant low- 
temperature and ground state correlations. Clearly, a 
second goal is to determine how to optimize particu- 
lar correlations in order to gain further insight into the 
mechanism responsible for the correlations, as well as to 
provide information which could be useful in the search 
for new materials with these correlations. Recently, the 
two-leg Hubbard ladder has been shown to exhibit power 
law d x 2_ y 2-\ike pair-field correlationsETH. Here we present 
a more detailed study of the dependence of the these 
pairing correlations upon the basic parameters of the 
model. We also examine their relationship to the low- 
energy single-particle spectral weight and the local anti- 
ferromagnetic correlations. 

The Hamiltonian for the two-leg Hubbard model with 
near-neighbor hopping is 



H 



\<t \ C i\a C i+l\a 



h.C.j -i_LEi, CT ( C L C i2<r 
+ h.C.) + C/E i\ n i\] n i\l- (1) 



Here c], creates an electron with spin a on the i th rung 
of the A th leg, with i = 1, . . . , L and A = 1 or 2. The 
intra-leg one-electron near-neighbor hopping matrix ele- 
ment is t and the inter- leg hopping is t±. The on-site 
Coulomb interaction is U. We will measure energies in 
units of t so that the basic parameters of the two-leg sys- 
tem become the hopping anisotropy t±/t, the ratio U/t of 
the interaction strength to the intra-leg hopping and the 
average electron filling (n) = l/(2£) Y^i,\( n i^ + n i*l)- 
At half-filling, with U/t > 0, the two-leg Hubbard lad- 



der is found to have both a charge gap and a spin gap 
Thus it is insulating with exponentially decaying short- 
range antiferromagnetic correlations. When the ladder is 
doped away from half- filling, the charge gap disappears, 
but a reduced spin gap remains for a range of doping and 
t±/t values. In addition, d x 2_ y 2-Hke power-law correla- 
tions are observedErti The internal structure of the pairs 
is illustrated in Fig. [j], in which we show the amplitude 



-rT C r'l C rl C r'T 



(2) 



for adding a singlet pair on near-neighbor sites along and 
across the legs, calculated using the density matrix renor- 
malization group (DMRG) method. Note the d x 2_ y 2-\ike 
change in sign of this matrix element. Here \Ni) is the 
ground state with four holes relative to the half-filled 
band and | JV2) is the ground state with two holes on a 
2x 16 ladder. For U/t = 8 and tj_/t = 1.5, the d x 2_ y 2-\ike 
structure of the matrix element, Eq. (||), extends over a 
region |r — r'| of order 4 rungs. 
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FIG. 1. Schematic drawing of the pair- wave function 
showing the values of the off-diagonal matrix element 
(-^2 1 (c^ j c^, . — c^c^^IA^i) for creating a singlet pair between 
near-neighbor sites. 



DMRG techniques^ have also been used to calculate 
the ground state expectation value of the rung-rung pair- 
field correlation function 



D(i,i) = (A(i)At(j)>. 



Here 



(3) 



(4) 



creates a singlet pair across the i th rung. For techni- 
cal reasons, the DMRG calculations are carried out for 
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2 x L ladders with open boundary conditions. For this 
reason, to obtain information on the spatial decay of the 
pairing correlations, we have averaged the pair-field cor- 
relation function D(i,j) over (typically) six pairs 
with I = \% — j\ fixed. This averaging procedure starts 
with symmetrically placed values and then proceeds 
to shift these to the left and right of center. When — 
approaches the lattice size, the number of possible 
pairs is reduced due to the proximity to the boundaries. 
Results for D{£ — computed in this way for 2 x 16, 

2 x 24 and 2 x 32 lattices with t x /t = 1.4, U/t = 8, and 
(ri) = 0.875 are shown in Fig. 0. The dotted line corre- 
sponds to a power law decay I . As can be seen, there 
are clear finite size effects when I approaches L. How- 
ever, one can determine when this occurs by comparing 
ladders of increasing length. In the following, we will 
discuss results obtained for 2 x 32 ladders over distances 
£ < 20, for which the end effects are negligible. 

In Fig. H we show D(£) versus £ for various values 
of tj_/t with U/t = 8 and (n) = 0.875. As this log-log 
plot shows, the pair-field correlations exhibit behavior 
consistent with a power-law decay 



0.1 



0.01 



0.001 



0.0005 b 



III 

1 


1 1 
• 


t x /t=l.o : 




o 


t x /t=1.2 " 




□ 


t x /t=1.3 - 




■ 


t x /t=l,4 _ 




▲ 


t x /t=1.6 


\ \c 






\ >v \ 






\ >»\ 
\ *A 






\ 11 
\ J 






2x32 






U/t=8 






(n) = 0.875 


\ 

\ 

\ 

\ 




— Ill 


l l 


\ aSml - 



/ 



10 



20 



(5) 



FIG. 3. D (£) versus £ for various values of t±/t with 
U/t = 8 and (nj = 0.875. 



for t±/t < 1.6. For t±/t = 1.6, D(£) shows large oscilla- 
tions and markedly reduced strength which we believe is 
associated with a transition to a phase in which the an- 
tibonding band is unoccupied. In Fig. |]we show the ex- 
ponent 9 versus t± /t obtained from a linear least-squares 
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FIG. 2. The rung-rung pair-field correlation function D(£) 
versus £ for U/t = 8, t±/t = 1.4, (n) = 0.875 on 2 x L ladders 
with L — 16, 24 and 32. The dashed and the dotted lines 
vary as £~ 2 and £~ 1 , respectively. 
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FIG. 4. Exponent 9 versus t±/t for U/t — 8 at fillings 
(n) = 0.75, 0.875 and 0.9375. 

fit of log D versus log(^) for I = 1 to 18 at various fill- 
ings. Because of the presence of some oscillations in D(£) 
and the finite range of £ available, the values of 9 thus 
obtained can only be rough estimates, but they do allow 
one to examine the trend in 9 as t±/t is varied. We find 
that 9 initially decreases with increasing t±/t until t±/t 
becomes of order 1.3 to 1.4 for U/t = 8 and the fillings 
shown. For larger values of t±/t, we find that the anti- 
bonding band is no longer occupied and the pairing cor- 
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relations rapidly collapse. Thus the pairing correlations 
are enhanced near the point at wh«jij|-|the antibonding 
band moves through the Fermi levefl'tffi. 

Another measure of the strength of the pair field cor- 
relations, which we will use, is the average of D(£)/D(l) 
for rung separations I = 8 to 12: 



D 



12 



(6) 



At a distance I — 8, we have moved beyond the cor- 
relation length characterizing the size of a pair and are 
probing the pair center-of-mass correlations. Figure ^ 
shows D versus t±/t for U = 8t at fillings (n) — 0.75, 
0.875, and 0.9375. Again, the pairing response is seen to 
initially increase with t±/t, reach a peak value, and then 
rapidly fall off. The variation of D with U/t is shown in 
Fig. Here we see that the peak value of D increases 
with U/t reaching a broad maximum for U/t of order 3 
to 8. In addition, the value of t±/t at which the peak oc- 
curs shifts towards smaller values as U/t increases. Also 
shown is D calculated for U — on the infinite ladder, 
showing the strong enhancement of the pairing correla- 
tions for the interacting system over those of the nonin- 
teracting system. In Fig. [t] we show results on how the 
minimum value of the exponent, # m i n , varies as a func- 
tion of U/t at (n) = 0.9375. Here, # min is obtained by 
varying t±/t while keeping U/t fixed. 

From the results shown in Figs. [| through 0, it is clear 
that the strength of the pairing correlations as measured 
by .D, Eq. (||), depend sensitively on t±/t and the filling 
(n) as well as on U / 1. The fact that D decreases for large 
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FIG. 5. D versus t±/t for U/t = 8 at fillings (n) 
0.875 and 0.9375. 
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FIG. 6. D versus t±/t for various values of U/t at filling 
(n) = 0.9375. 
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FIG. 7. Minimum value of the exponent, (9 m i n , versus U/t 
at filling (n) = 0.9375. 

values of U/t is consistent with the variation of the ex- 
change interactions J = At 2 /U and J± = itj_/U, which 
become weaker at large values of U. It is interesting to 
note that D has its largest value for intermediate coupling 
U /t. Since J±/J varies as (t±/t) 2 , the interchain antifer- 
romagnetic correlations are enhanced relative to the in- 
trachain correlations as t±/t increases. The z—z spin cor- 
relation function (M^MJ A ), where M* x = Jij,At ~ n i,\U 
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is shown for interchain (J = i, A = 2) and intrachain 
(j = i + l, A = 1) nearest neighbor sites in Fig. ||. Thus 
both U/t and t±/t enter in determining the strength and 
anisotropy of the exchange couplings and the local struc- 
ture of the antiferromagnetic correlations on the ladder. 
The underlying pairing interaction is clearly associated 
with these short-range antiferromagnetic correlations. 
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FIG. 8. The inter- (j = i, A = 2) and intra-chain 
(j = i + 1, A = 1) near- neighbor magnetic correlation function 



tj_/t for (n) = 0.875 and = 8. 



In addition, and of particular importance, the ratio 
t_i_/i) the filling (n), and the interaction strength U/t 
determine the energy and momentum structure of the 
single-particle spectral weight 



A(k,u) = -~ImG(k,uj) 

7T 



(7) 



Results for A(k, u), obtained from a maximum entropy 
analytic continuation of Monte Carlo data for a 2 x 16 
lattice with periodic boundary conditions, are shown in 
Fig. | for U/t = 2 and in Fig. 10 for U/t = 4.0 Here, 
we have chosen (n) = 0.94, the filling at which the pair- 
ing correlations have maximum strength in Fig. ||, and 
a temperature T = 0.125i, which is sufficiently low to 
ensure that the band structure would not shift at still 
lower temperatures. (Note that (n) is a measured quan- 
tity in the grand canonical Monte Carlo simulations, and 
is thus not fixed as in the DMRG calculations for which 
the number of particles is definite.) We show data for the 
k± = (bonding) and k± = 7r (antibonding) branches su- 
perimposed in a density plot in which the density of the 
shading represents the relative amount of spectral weight 
and is plotted as a function of the momentum along the 
chains, k, and the energy u>. Although the resolution is 



limited by the finite temperature and statistical errors, 
one can see coherent, dispersive bonding and antibonding 
bands. 
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FIG. 9. Single-particle spectral weight A(k,uj) versus ui at 
different values of k for U/t — 2, (n) = 0.94 and a tempera- 
ture T/t = 0.125. Here results are shown for (a) tx/t = 1.6, 
(b) 1.8 and (c) 2.0. The bonding (k± = 0) and antibonding 
(k± = 71") branches are superimposed on the density plot, with 
the indicated symbols marking the positions of peaks associ- 
ated with each band. Areas of darker shading correspond to 
areas of higher relative spectral weight. The dotted lines in- 
dicate the position of the noninteracting (U = 0) bands and 
the solid line at ui = indicates the Fermi level. 
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The spectral weight distribution clearly evolves with 
t±/t for both values of U/t. For U/t = 2, the width and 
dispersion of the bands follow quite closely those of the 
nonintcracting, [7 = 0, bands, which are indicated by 
dotted lines in Fig. [| As tj_/t increases from 1.6 to 2.0, 
the peaks in the bonding and antibonding spectral weight 
at the Fermi surface move towards a momentum separa- 
tion Ak = (-7T, tt) where scattering from the short range 
antiferromagnetic correlations can become most effective. 
The point at which the bottom of the antibonding band 
is just at the Fermi level, t±/t = 1.8, (Fig. ||(b)) agrees 
well with the position of the peak in D in the U/t = 2 
curve in Fig. [| At tyjt = 1.8, there is large amount of 
single-particle spectral weight near the Fermi level in the 
bonding band near k = tt and in the antibonding band 
near k = 0. At larger values of t±/t, the antibonding 
band pulls away from the Fermi energy, leaving only the 
bonding band with spectral weight at the Fermi energy, 
as shown in Fig. |9|(c) for t±/t = 2.0. We believe that 
it is the variation in spectral weight with t±, coupled to 
antiferromagnetic fluctuations which are strongly peaked 
near (71", tt) that is primarily responsible for the peak in 
D versus t±/t 

As shown in Fig. |^, as the on-site coupling is increased 
to U/t — 4, the value of t±/t at which D peaks shifts 
to smaller values. This is due to narrowing of the quasi- 
particle bands relative to the U = bands for larger U. 
In Fig. [l(], we exhibit the single-particle spectral weight 
at t±/t = 1.4, 1.6, and 1.8, values which bracket the peak 
in D, which occurs at t±/t = 1.6. As in Fig. |j| one can 
clearly see that the maximal spectral weight at the Fermi 
level in the antibonding band at k = and the bonding 
band at k — tt occurs just at the value of t±/t associated 
with the peak. One can also see additional structures 
in both bands, which are reflections of the portion of 
each band below the Fermi level about the Fermi level. 
These structures, present due to stronger antiferromag- 
netic correlations at the larger U ft value, are remnants 
of features that would be generated by the halving of the 
Brillouin zone in an antiferromagnetically ordered state. 
Both bands are flattened relative to the U = bands 
near the Fermi level, i.e. at k = in the antibonding 
band and at k = tt in the bonding band, especially for 
t±/t = 1.6 which is shown in Fig. |l^(b). As a result 
of this, there is large amount of single-particle spectral 
weight near the fermi level. This behavior is similar to 
the buildup of spectral weight near (tt, 0) observed in sim- 
ulations of the n two-dimensional Hubbard model doped 
near half-fillingcffl. 

These results show that near half-filling, the strength 
of the pairing correlations depends sensitively upon tj_ jt 
and the band filling as well as U/t. The optimum value of 
D occurs for intermediate values of U / 1 and for (n) near 
half- filling with t±/t ~ 1.5. We find that for these val- 
ues, there are strong inter-chain antiferromagnetic cor- 
relations and a large single-particle spectral weight at 
the Fermi surface points of the bonding and antibonding 
bands. 
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FIG. 10. Single-particle spectral weight A(k, uo) versus k 
and uj plotted as in Fig. 9 for the same parameters as in Fig. 
9, except with U/t = 4 and (a) tx/t = 1.4, (b) 1.6 and (c) 1.8. 
The dotted lines indicate the position of the noninteracting 
(U = 0) bands and the solid line at uj — indicates the Fermi 
level. 
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